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' We overview the general scaling behavior of the effective temperature in the quenches of simple 

, non disordered systems, like ferromagnets, to and below Tc- Emphasis is on the behavior as di- 

mensionality is varied. Consequences on the shape of the asymptotic parametric representation are 
^ ' derived. In particular, this is always trivial in the critical quenches with Tc > 0. We clarify that the 

Q quench to Tc = at the lower critical dimensionality d_t, cannot be regarded as a critical quench. 
Implications for the behavior of the exponent a of the aging response function in the quenches below 
\^ I Tc are developed. 
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^ ■ I. INTRODUCTION 

C/5 , 

In the study of slow relaxation phenomena [1], typically, a temperature quench is made at time t = and the 
quantities of interest, such as the autocorrelation C{t,tw) and the linear response function R{t,tyj), are monitored 
at subsequent times Q < < t. The peculiar and interesting feature of these processes is that equilibrium is never 
' reached, namely there does not exist a finite time scale teq such that for > teq two time quantites become time 
^ I translation invariant. 

O ' A substantial progress in the study of these phenomena has been achieved after Cugliandolo and Kurchan [2] have 
, ^ , introduced the modification of the fluctuation-dissipation theorem (FDT) as a quantitative indicator of the deviation 
from equilibrium [3]. This can be formulated in several equivalent ways. Here, we find convenient to use the effective 
T— I , temperature [4] 

> ■ 

(N 



which coincides with the temperature T of the thermal bath if the FDT holds, that is in equilibrium, while it is 
different from T off equilibrium. Because of the absence of a finite equilibration time, the characterization of the 
' system staying out of equilibrium for an arbitrarily long time is contained in Imit^^ooTef f{t,tw)- However, since 
"Y"^ tit, — !■ oo implies also f — s- cx), to make the limit operation meaningful one must specify how the two times are pushed to 
infinity. One way of doing this is to fix the value of the autocorrelation function C(i, t^) as oo. More precisely, 

since C(t, t^) for a given tw is a monotonously decreasing function of t, the time t can be reparametrized in terms of 
^ C obtaining T^f f(t,tw) = Tf,f f{C,tw) and, keeping C fixed, the limit 

O : r(C)= hm fe//(C,iJ (2) 

o 

defines the effective temperature in the time sector corresponding to the chosen value of C . This quantity is important 
because, when appropriate conditions are satisfied, provides the connection between dynamic and static properties [7] 
through the relation 

P{q) = T-^T{Cr%^q (3) 

where P{q) is the overlap probability function of the equilibrium state. Therefore, different shapes of T(C) are 
associated to different degrees of complexity of the equilibrium state. A first broad classification of systems has been 
made [8,9] drawing the distinction between coarsening systems, structural glasses and the infinite range spin glass 
model on the basis of the patterns displayed by T{C), or related quantities such as the fluctuation-dissipation ratio 
(FDR) and the zero field cooled (ZFC) susceptibility. In this paper we focus on coarsening systems, exemplified by 
a non disordered and non frustrated system such as a ferromagnet, relaxing via domain growth after a quench to 
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or below the critical point [10-12]. Wc show that even within the framework of these so called trivial systems, the 
spectrum of the behaviors of T{C) is quite rich and interesting. 

In order to illustrate the problem, in Fig.l we have drawn schematically the phase diagram, in the temperature- 
dimensionality plane. The critical temperature vanishes at the lower critical dimensionality dh- The disordered and 
ordered phase are, respectively, above and below the critical line Tc{d). Slow relaxation arises for quenches with the 
final temperature everywhere in the shaded area, including the boundary. Usually, one considers some fixed d > 
(dashed line in Fig.l) and studies the quench to a final temperature T <Tc- The form of T(C) in the critical quench 
with T = Tc (Fig. 2a) is different from that with T < Tq (Fig. 2b) and both are independent of d. The interesting 
question, then, is what happens at the special point {d^^T = 0), where critical and subcritical quenches merge. 
Namely, will T(C) display continuity with the critical or with the subcritical shape? On the basis of the available 
evidence from analytical solutions [13-16] and numerical simulations [17] , the surprising answer is with neither of 
them. A third and nontrivial form of T(C) is found (Fig. 2c) in the quench to (d^, T = 0). 

In this paper we present the scaling framework which allows to account in a unified way for this wealth of behaviors 
displayed by T(C) upon letting T and d to vary. This, in turn, allows to gain insight on the challanging problem [16,18] 
of the scaling of R{t, tw) in the quenches below Tc- Furthermore, on the basis of these ideas, predictions can be made 
for the interesting case of the XY model. In that case, the phase diagram of Fig.l is enriched by the presence of 
the Kosterlitz-Thouless (KT) line of critical points at the lower critical dimensionality (d/, = 2), with T < Tkt- 
Therefore, slow relaxation in the XY model arises also by arranging the quenches onto the KT line. Since this is a 
critical line, the phenomenology for < T < Tkt is expected to be akin to the one along the Tc{d) line, while at 
T = the switch to the third nontrivial form of T(C) is expected to take place. 

The paper is organised as follows. In section II the scaling properties of C{t,tyj) and R{t,t^) are reviewed. In 
section III the general scaling behavior of Tgff{t,tyj) is derived, and in section IV the parametric form T(C) is 
obtained. The problem of the exponent a of R{t, t^) below Tc is presented in section V. Sections VI and VII 
are devoted to the illustration of the general concepts in the context of the large- iV model and of the XY model, 
respectively. Conclusions are presented in section VIII. 

II. SCALINGS OF C{T,Tw) AND R{T,Tw) 

In the following we consider quenches from an initial disordered state above the critical line to any one of the states 
in the shaded area of Fig.l. For definiteness, we shall refer to purely relaxational dynamics without conservation of 
the order parameter, but the results are of general validity. The common feature of all these processes is that, after 
a short transient, dynamical scaling sets in with a characteristic length growing with the power law L{t) ~ t^^^ [10]. 
The scaling properties with T = Tc are different from those with T < Tc [11]. In the renormalization group language 
this means that for a given d there are two fixed points, one unstable at Tc and the other stable and attractive at 
T = (thermal fluctuations are irrelevant below Tc) [19]. The dynamical exponent z on the Tc{d) line coincides with 
the dynamical critical exponent Zc and depends on dimensionality, while below criticality it takes the dimensionality 
independent value z = 2. These values become identical at {dL,T = 0), since lim^j^di Zc = 2 [20]. 

Let us, next, review the scaling properties of C{t, t^) and R{t, t^)- The fixed point structure enters in two ways: in 
the values of the exponents involved and, less obviously, in the very same form of the functions C{f., t^,) and R{t, tw)- 
In analysing the asymptotic behavior of these quantities, it is necessary to distinguish between the short time and the 
large time behaviors obtained by letting the waiting time t^ to become large, while keeping, respectively, either the 
time difference t = t — tw ov the time ratio x = t/tw > 1 fixed. Notice that from x = 1 — r/tw follows that the short 
time regime, when using the x variable, gets compressed into x = 1. 

An important general property of two time quantities in slow relaxation phenomena is that in the short time regime 
the system appears equilibrated and the FDT is satisfied, while the off equilibrium character of the dynamics, or 
aging, shows up in the large time regime [1,3,12]. The mechanism underlying this property, however, is quite different 
for quenches to and below Tc- 

A. Critical quench: multiplicative structure 

In the critical quenches, for tw sufficiently large, C{t,ty,,to) and R{t,tyj,to) take the product forms [11,12] 

C{t, tw, to) = (r + to)-''gc{x, y) (4) 

R{t, tu,, to) = {t + to)-^'+''^gR{x, y) (5) 
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where gc,R{x, y) are smooth functions of x with a weak dependence ony = ta/tw In addition to the observation times, 
t and tw, we have exphcitely included also the dependence on a microscopic time to, which is needed to regularize 
these functions at equal times. In particular, we shall take C{t,t,to) = 1 throughout. In the short time regime 

C{t,t^,to) = {T + to)-''gc{l,0) (6) 

R{t,UM) = (r + io)-('+"^ffK(l,0) (7) 

which are the autocorrelation and response function of the stationary critical dynamics satisfying, therefore, the 
equilibrium FDT, i.e. TcR{t) = —drC{T), which implies 

TcgR{l,0) = bgc{l,0) (8) 

and 

a = b. (9) 

We emphasize that, in the critical quench, as a consequence of the multiplicative structure and of the FDT, the 
exponents a and b are not independent and coincide. Furthermore, their common value is given by [11,12] 

a = b= {d-2 + 'n)/zc = 2()/uZc (10) 

where r?, (3 and v are the usual static exponents. Using the geometrical interpretation of the critical properties Eq. (10) 
can be rewritten as a = = 2{d — D)/zc, where D is the fractal dimensionality of the correlated critical clusters [21]. 
These become compact as Tc ^ yielding 

lim a = 6 = 0. (11) 

On the KT line, where d = 2, Eq. (10) is replaced by a = 6 = 'n{T)/z with r?(T) vanishing as T ^ [22]. Eqs. (4) 
and (5) can be recast in the simple aging form 

C{t,t^,to)=t-^fc{x,y) (12) 

R{t,t^,to)=t-^^+''^fR{x,y) (13) 

where 

fc{x, y) = {x-l + y)~^gc{x, y) (14) 

fR{x, y) = {x-l+ y)-(i+»)5fl(a;, y). (15) 
These functions, for large x and tw, decay with the same power law 

fcAx,^) = Ac,RX-^^''^ (16) 
where Ac is the autocorrelation exponent [11,12]. 



B. Quench below the critical line: additive structure 

Below Tc, the picture is different because the stationary and aging components enter additively into the autocor- 
relation function [23,1,3] 

C{t,tu,,to) = Cstir) + Cag{t,tu,,to) (17) 

where Cst(T) is the equilibrium autocorrelation function in the broken symmetry pure state at the temperature T. 
Prom the equal time properties 

C{t,t,to) = l Cstir = 0) = 1 - (18) 
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which imply 

Cag{t,t,to) = (19) 

where M is the spontaneous magnetization, follows that in Eq. (17) the stationary component corresponds to the quasi- 

equilibrium decay to the plateau at the Edwards- Anderson order parameter qea = M"^, while the aging component 
Cag describes the off-equilibrium decay from the plateau at much larger times. In the short time regime, then, we 
have 

CiT + t^,t^,to) = Cst{T)+M'^ (20) 
while, taking into account that C'st{T) vanishes as r — > oo, in the large time regime 

C{t,ty),to) = Cag{t,tu],to)- (21) 

From this it is easy to verify that, contrary to what happens with the multiplicative structure of Eq. (4), the limits 
tw ^ oo and f — > oo do not commute 

lim lim C{t,t^,to) = M"^ (22) 



lim lim. C{t,t^,to) =Q (23) 

yielding the weak crgodicity breaking [1], characteristic of the quenches below Tc- 
The additive structure of the response function 

R{t, t^, to) = Rstir) + Rag{t, t^,to) (24) 

is obtained in the following way [3]: in the short time regime stationarity holds R{t,tyj,to) = Rst{T) and Rst{T) is 
defined from the FDT 

Rstir) = -^§^Cst{T) (25) 

whilst Rag{t,tyj,to) remains defined by Eq. (24). This implies that Rag{t,tw,to) vanishes in the short time regime 
and since, conversely, Rst{T) vanishes in the large time regime, we may write 

R{t,t^,to) = I R^^^llt^^t^' for" \>1. (26) 

The components Cag{t,twTtQ) and Rag{t,t-uj,tQ) obey simple aging forms like (12) and (13) with scaling functions 
fc,R{x,y) decaying, for large x, with a power law of the form (16). It is understood that all the exponents must be 
replaced by their values below Tc, which are different from those at Tc- In particular, from phase ordering theory it 
is well known [10] that 6 = independently from dimensionality. This can be readily imderstood from h = 2{d — D)/ z, 
which holds in general, and from the compact nature of the coarsening domains, which gives D = d. Quite different 
is the case of the exponent a. Since the FDT relates only the stationary components Cst{T) and Rstir), due to the 
additive structure of Eqs. (17) and (24) we have no more the constraint a = 6, as in the critical quench. Namely, in 
the quench below Tc the value of a is decoupled from that of b. As a matter of fact, the determination of a below Tc 
is a difficult and challanging problem [17,18], which will be discussed in section V. 



C. Quench at (c?z,,T = 0): additive structure 

As we have seen above, the structures of C(t, tw,to) and i?(t, tw,to) follow different patterns in the quenches to and 
below Tc- This raises the question of which one applies when the quench is made at (d^.T = 0). In other words, 
should we consider the state idL,T = 0) as a critical point with Tc = 0, as it is often done in the literature, or should 
we consider it as the continuation to rfj, of the T = states below criticality? The correct answer is the second one, 
since at d^ the system orders, as in all the states below the critical line, while on the 7c (d) line there is no ordering. 
In particular, this implies that weak ergodicity breaking takes place also in the quench to (d/,, T = 0). From Eqs. (17) 
and (18) then follows Cst(r) = and 
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C{t,tyj,to) = fc{x,y) 



(27) 



with /c(l, 0) = 1, since = 1 at T = 0. 

The fact that the quench to {dL,T = 0) belongs to the additive scheme becomes important when considering the 
response function, because if Cstir) vanishes in the ground state, the same is not true for Rgtir). Therefore, Eq. (24) 
holds with Rst{T) ^ and it is important to subtract this contribution from R{t,ty^,,to) if one wants to study the 
scaling properties of Rag{t, tw, to)- We shall come back on this point in section VI, when treating the explicit example 
of the large- A?' model. It should be mentioned, here, that the linear response function at T = is well defined only 
for systems with soft spins. In the case of hard spins there docs not exist a linear response regime when T = 0. For 
a discussion of this problem and how to bypass it in Ising systems, we refer to [13,16]. 



D. Comparison with theory 

A natural question is how the picture outlined above compares with theoretical approaches. In the case of the 
quench to Tc ronormalization group methods are available [24,12], which account for the multiplicative structure 
of the autocorrelation and response function and give explicit expressions for the quantities of interst up to two 
loops order. The development of systematic expansion methods for the quench below Tc is much more difficult [25]. 
Recently, Henkel and collaborators [26,27] have developed a method based on the requirement of local scale invariance. 
According to this approach the response function obeys the multiplicative structure (5) both in the quenches to Tc 
and below Tc- Problems, then, arise in the latter case. The main one is of conceptual nature, since, as explained 
above, the weak ergodicity breaking scenario requires the autocorrelation function to obey the additive structure (17). 
Then, in the short time regime, where Eq. (20) holds, from the FDT one has dt^Cst{T) ^ {t + tQ)^'^^^°-\ which implies 
Cstij) ~ (r + to)""- In other words, if the local scale invariance prediction for the response function is valid then the 
stationary correlation function must decay with a power law. This is the case, as we have seen above, in the critical 
quenches and, as we shall sec in section VI, in the large A'^ model. Conversely, in all other cases where the equilibrium 
correlation function below Tc decays exponentially, the response function cannot be of the form (5). Nonetheless, 
Henkel et al. do criticize [27] the additive structure below Tc- 

The second problem is that existing results, exact [13-15] and approximated [16,28,25] as well as numerical [18,29], 
lead to an the aging part of the response function which, for systems with scalar order parameter [30] , is of the form 

Rag{t, t^, to) = t-'"'{T + to)-^''-^''+^^gR{x, y)- (28) 

This is qualitatively different from the multiplicative form (5), both in the structure and in the exponents. As 
mentioned above, an exception is the large N model, where R{t,tw,to) is of the form (5) also below Tc- Yet, as we 
show in section VI, the additive structure applies to the large N model as in all cases of quench below Tc- 

III. EFFECTIVE TEMPERATURE 

The next step is to see how these different behaviors of C{t,tu],to) and R{t,tw,to) do affect the behavior of 

Teff{t,tu,,to)- 

A. Critical quench 

Prom Eqs. (12), (13) and (9), using the definition (1), we get 

Teff{t,t^,to)=F{x,y) (29) 

with 

F{x, y) = -fecix, y)/fR{x, y) (30) 

and 



fdc{x,y) = bfcix,y) + 



x-^fc{x,y)+y-^fcix,y) 



(31) 
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In the short time regime, Eqs. (6), (7) and (8) lead to the boundary condition 

F{l,y)=Tc 

while, in the aging regime, from Eq. (16) follows 
where Xoo is the limit FDR of Godreche and Luck [11]. 



(32) 



(33) 



B. Quench below the critical line 

Below Tc, where b = 0, from the definition (1), the additive structure and Eq. (25) we get 



Teff{t,tw,to)) 

which can be rewritten as 



Rst{T) + Rag{t,tw,to) Rst{T) + Rag{t,tw,to) 



Teff{t,tw,tQ)) = T 



+ tlF{x,y) 



Rag (^1 iw , ^o) 



Rst{T) + Rag{t,tui,to) 



(34) 



(35) 



where F{x,y) is still defined by Eqs. (30), (31) with fc{x,y), fR{x,y) the scaling functions of Cag{t,tyj,to) and 

Rag (^1 , ^o) • 

Using Eq. (26) we may approximate Eq. (35) with 



where 



nff{t,t^,to)) = TH{x)+tlF{x,y)[l - H{x)] 



(36) 



^(-) = |0 for :>1. (37) 
Therefore, for large tyj the general formula containing both critical and subcritical quenches is given by [31] 

( T for a: = 1 

Teff{t,tw,to)) = ^^a^-l,p^^^y^ for X>1. ^^^^ 

Of course, this includes also the quench to {cIl, T = 0). 



IV. PARAMETRIC REPRESENTATION 

In order to derive the parametric representation, as explained in the Introduction, we must express the time 

dependence of Teff{t,tw,tQ) through C{t,tw,tQ) and then let tw oo, obtaining T{C). Prom Eqs. (4) and (17), for 
large , the general form of the autocorrelation function can be written as 



C{t, twjto) 



1 for X = 1 

tw^fc{x,y) for x>l. 



(39) 



The task is to eliminate x between Eqs. (38) and (39), which now we do separately for the quenches into the different 
regions of the phase diagram. 
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A. Critical quench 



In this case a = b > 0. Prom Eqs. (12) and (14), letting tw ^ oo with x fixed, we obtain the singular limit 

whose inverse is readily obtained exchanging the horizontal with the vertical axis 

f oo for C = 
^(^) = |l for 0<C7<1. (41) 

Inserting into Eq. (38) we obtain (Fig. 2a) 

^(^)-\Tc for 0<C<1 (42) 

where we have used Eqs. (32) and (33). This is a universal result, since all the non universal features of F{x, 0) for the 
intermediate vahics of .t have been cUminatcd in the limit process. Therefore, wc have that for all quenches to Tc > 0, 
apart for the value Too at C = 0, the parametric plot of T(C) is trivial in the sense that the effective temperature 
coincides with the temperature Tc of the thermal bath. Notice that this implies that for the ZFC susceptibility 

X{t,t^) = j^ dsR{t,s) (43) 

whose parametric form is related to T(C) by 

'=(^) = XrCT (") 

a linear plot is obtained 

X{C) = (1 - C)/Tc (45) 

as for equilibrated systems [3]. 



B. Quench below the critical line 

In this case 6 = and from Eq. (39) follows 

^(^) = { /c(x,oT fr^>i (46) 

where, recalling Eqs. (18) and (21), fc{x,0) is the smooth function describing the fall below the plateau at M^. 
Therefore, the inverse function is given by 

^^^^~\1 for M2<C<1 

and inserting into Eq. (38) we find 

f t )-( tinfc\C),0) for C<M^ 
Je//lO,t^j ^ j y for AP<C<1. 

Although the actual value of a, as will be explained in section V, is to some extent a debated issue, there is consensus 
on a > for d > (11- Therefore, taking the tyj ^ oo limit we recover the well known result [1,3] (Fig. 2b) 
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Again, all the details of F{x, 0) having disappeared, the function T(C) is universal. 

Let us remark, here, as a further manifestation of the difference between quenches to and below Tc, that the 
behavior (42) of T{C) cannot be recovered by letting T — > Tq and M — > in (49). In fact, in the latter case we find 

^limT(C;) = {- (50) 
which differs from Eq. (42), where Too is a finite number. 

C. Quench to (dt,T = 0) 
If we take the limit d ^ dz, in Eqs. (42) and (49) we obtain, respectively 

for 0<C<1 (^1) 

and 

which are very much different one from the other. So, there is the problem of which is the form of T(C) in the 
quench to (di, T = 0). A statement to this regard can be made on the basis of the substantial amount of information 
by now accumulated from sources as diverse as the exact solutions of the Id Ising model [13,14] and of the 2d 
large- model [15], the Ohta-Jasnow-Kawasaki type of approximations with d = 1 [16], in addition to numerical 
simulations [17] of systems at rfj, with scalar and vector order parameter, with and without conserved dynamics. 
In all of these cases we have obtained the parametric plot (44) of the ZFC susceptibility, which gives for T(C) at 
(di, T = 0) a non trivial and non universal finite function of the type depicted in Fig. 2c. By contrast, there does not 
exists, up to now, any evidence for behaviors of T{C) of the type (51) or (52). 

Assuming that the generic behavior of T(C) is of the type in Fig. 2c, as the amount of evidence quoted above 
strongly suggests, this is compatible with Eq. (38) only if a = 0. Then, from Eqs. (46), (47) and (48) with = 1 
we may write 

T{C)=F{fc\C),Q) (53) 

which represents a smooth, finite, non trivial function deacreasing smoothly from Too = -F(oo,0) at C = toward 
zero at C = 1, preserving all the non universal features of fcix, 0) and F{x, 0) for the intermediate values of C. 

D. Physical temperature and connection with statics 

We conclude this section with comments on the interpretation of T(C) as a true temperature and on its connection 
with the equilibrium properties of the system. 

The relation of the effective temperature with a physically measurable temperature is a very interesting issue, 
which has been thoroughly investigated in Refs. [4-6]. A necessary condition for the identification of the effective 
temperature with a true temperature is its uniqueness, that is the indcpence from the observable used in the study of 
the deviation from FDT. For what concerns the quenches to Tc, in [6] it is shown that the multiplicative structure (4) 
and (5) is observable independent. Therefore, the derivation of Eq. (42) is also observable independent, except for the 
actual value of T^a-, which does to depend on the observable [6]. In this case, the criterion of uniqueness does indeed 
lead to the identification of the effective temperature with a physical temperature, since T{C), for C > 0, coincides 
with the temperature of the thermal bath Tq- Conversely Too, due to the non uniqueness of its value, cannot be 
considered as physical temperature. Although results of comparable generality are lacking for the quenches below Tc, 
we expect also Eq. (49) to be observable independent, due to the universality of T(C) (for observables with a > 0). 
Finally, the effective temperature (53), in the quench to {dL,T = 0), it is certainly not physical, since, as explained 
above, it does to depend on the details of the scaling function, which makes it observable dependent. To this case 
belong the results for the d = 1 Ising model [5], which do indeed show the observable dependence of the effective 
temperature. 
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The connection between static and dynamic properties is given by Eq. (3). We can now assess the status of 
this relation, in the hght of the results derived above. In all cases, the overlap probability function is given by 
P{q) = S{q — AP). We must, then, verify if the r.h.s. of Eq. (3) actually is a (5 function. Prom Eqs. (42,49,53) 
follows that T-^T{C)~^\c=q vanishes everywhere, except at C = M^. Therefore, in order to establish whether this 

is a (5 function, we must look at the integral dC tJc) ' This gives finite values in the quenches to Tc and below 
Tc- The validity of Eq. (3), instead, cannot be established in the quench to {dL,T = 0). In that case the value of 
the integral is not determined, since the contribution at the upper limit of integration is given by the ratio of two 
vanishing quantities. 

V. THE RESPONSE FUNCTION EXPONENT 

With the results for T{C) illustrated in the previous section we come to grips with the problem of the exponent 
a. The vanishing of a, required by the behavior of T{C) in the quench to {d^, T = 0), cannot be accounted for by 
Eq. (11) because, as explained in section II C, the quench to {d^^T — 0) is not a critical quench. Hence, the vanishing 
of a must be framed within the behavior of a in the quenches below the critical line. 

According to a popular conjecture [32], a for T < Tc should coincide with the exponent n/z entering in the time 
dependence of the density of defects, which goes like with n = 1 or n = 2 for scalar or vector order parameter, 

respectively [10]. We recall that for T <Tc the dynamical exponent z = 2 does not depend on d. Thus, in the scalar 
case one ought to have a = 1/2 and in the vector case a = 1, independently from d and, therefore, also at di,. This is 
obviously incompatible with the large body of evidence for a = at d^, quoted in the previous section. 

This difficulty is supersed in the alternative proposal which we have made [16,18] for the behavior of a as d is varied, 
on the basis of the exact solution of the large- AT model [15] for arbitrary d and the Ohta-Jasnow-Kawasaki type of 
approximation [16,28,25], also for arbitrary d. In these two cases one finds that a does depend on d according to 

where djj is a parameter dependent on the universality class. We have proposed to promote this result to a phe- 
nomenological formula, whose general validity we have tested by undertaking a systematic numerical investigation of 
systems in different classes of universality and at different dimensionalities [16-18]. The results obtained are in quite 
good agreement with Eq. (54), taking du = 3 oi du = 4, respectively for scalar [33] or vector order parameter. We 
emphasize that Eq. (54) gives limd^di a = when the limit is taken below Tc, and this should not be confused with 
Eq. (11) holding for the critical quenches. 

The final remark is that, up to now Eq. (54) remains a phenomenological formula. The challange is to make a 
theory for it. A preliminary attempt to relate the dimensionality dependence of a to the roughening of interfaces in 
the scalar case has been made in [17]. However, the explanation of Eq. (54) as a result of general validity seems to 
require an understanding of the response function much deeper than we presently have. 

VI. LARGE-iV MODEL 

In this section all the concepts introduced above are explicitely illustrated through the exact solution of the large- A^ 
model. This model, or the equivalent spherical model, has been solved analytically in a number of papers [34-36]. 
Here, we follow our own solution in Ref. [15], where we have shown that for quenches to T < Tc the order parameter 
can be split into the sum of two stochastic processes 4>{x, t) = a{x, t) + il^{x, t) which, for t sufficiently large, are 
independent and mimic the separation between the interface and bulk fluctuations taking place in domain forming 
systems [37]. Correspondingly, the autocorrelation function decomposes into the sum of the autocorrelations of a and 

C{t, tyj, to) = Ccr{t, tyj, to) + C^{t, tyj, to) (55) 

which obey the scaling forms 

Ca{t,tyj,to) = t-^'U{x,y) (56) 

and 

C^{t,tyj,to) = t-''*U{x,y). (57) 
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For the purposes of the present paper it is sufficient to limit the analysis to d < 4, where the exponents and the 
scaling functions are given by 



d-2 for T = Tc 
for T < Tc 



b^ = {d- 2)/2 



{x + l + y) 



-d/2 



IT 



d/2 - 2 for T = Tc 
-d/2 for T <Tc 



(58) 
(59) 

(60) 
(61) 

(62) 



and A„ is a constant, which coincides with in the quenches below Tq. The scaling function ft),{x,y) can also be 
rewritten as 



U{x,y) = {x-l + y) ''*g^{x,y) 



with 



T 



2/{x-l+y) 



/2 / dzil+zy'^ 



which allows to recast C.^{t,tyj,to) in the multiplicative form 

C^it,t^,to) = {t + tQ)-^^g^{x,y). 

The response function is given by 

R{t,t^,to) =1-^^+""^ fR{x,y) 

with 

fnix, y) = {Ai^r'^'^x - 1 + y)-^^+-) x-^ 

and 

a = h-^, = (d-2)/2. 
Finally, in the phase diagram of Fig.l the critical line is a straight line 

Tc = {ATTY'^K^-'^{d-dL)/2 

where A is the momentum cutoff and dL = 2. 



^-\{x- ^ + y)z 



(63) 

(64) 

(65) 
(66) 
(67) 
(68) 
(69) 



A. Scalings of C{t,t.w,tQ) and R{t,tw,to) 

Let us now extract from the above results the scaling properties of C(t,tw,to) and R{t,tw,to) on the different 
regions of the phase diagram. 
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1. Critical quench 

The first observation is that the autocorrelation and response function have the multiplicative structures (4) 
and (5). For R{t,tw,to) this is immediately evident from Eq. (67), where fR{x,y) is of the form (15) with 
gR{x,y) = {4w)-''/^x'-^/\ 

For C{t,tw,to), from b^r = 26^ follows that for large tyj the first term in Eq. (55) is negligible with respect to 
the second one, yielding 

C{t,tyj,to) = Cjj,{t,tyj,to) (70) 

and this is of the form (4) and (12) with the identifications b = b^, fc{x,y) = fip{x,y) and gc{x,y) = g^{x,y). 
Furthermore, from Eqs. (59) and (68) follows 

a = b={d-2)/2 (71) 

in agreement with Eq. (10), since in the large- A'^ model rj = and Zc = 2. Finally, from Eq. (64) we have 
gp(l,0) = (47r)-''/22rc/(d- 2) and, using 5k(1,0) = (47r)-'^/^ it is easy to check that Eq. (8) is verified. 

2. Quench below the critical line 

For large tw, the additive structure of Eq. (17) is found ready made in Eq. (55), with the identifications 

Cstir) = C^(r) = {T + to)-''^g^{l,0) (72) 
where g^{l,0) = {'i-K)-'^/^2T/{d-2) and 

Cag{tityj,to) = Ccr{t,tw,to) (73) 

which gives b = b^ = and ,fc{x, y) = faix, y). The pawer law decay (72) of the stationary component Cst{T) is 
a peculiarity of the large- model, since below Tc to lowest order in 1/A'^ only the Goldstone modes contribute 
to thermal fluctuations [38]. 

Carrying out the prescription outlined in section II B for the construction of the corresponding components of 
the response function, from Eq. (25) we have 

R,t{T) = (47r)-'^/2(r + to)-^'+^*^ (74) 
and using the identity b^ = a we find 

Rag{t,tw,to) = R{t,tw,to) — Rst{T) (75) 

= t-^'+''^fR{x,y) 



with 



fR{x,y) = {An)-''/^ 



x'i/^ - 1 
{x-l + ?/)i+"' 



(76) 



3. Quench to T = with d = 

As argued in Section V, this case belongs to the additive scheme and is contained in the previous one. From 
Eq. (61) follows that C^{t,ty,,to), and therefore also Cgtir), vanish identically in the quenches to T = 0. Then 
we have C{t,tiu,to) = C„{t,tyj,tf)) which, using Eqs. (56) and (58), is equivalent to Eq. (27). Conversely, Rsi{t) 
and Rag{t,tw,to) are temperature independent and are still given by Eqs. (74) and (76) with d = 2 and a = 0. 
Notice that the equilibrium FDT is satisfied also in the limit T ^ 0, where Cgtir) vanishes, since from Eq. (72) 
follows 



We emphasize, here, that the explicit form (68) of the exponent a is a particular case of Eq. (54) with n = 2. z = 
2,dL = 2 and du = 4. It should also to be remarked that a is given by the same expression (68) both in the quenches 
to Tc and below Tc- This is a peculiarity of the large- A'' model, where ri = and Zc = z. 
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B. Effective temperature T(C) 



1. Critical quench 

Prom the explicit expressions for fc{x,y) and fR{x,y) in the critical quenches follows 

Teff{t,t^,to) = F{x,y) (78) 
= dTc {(^ + 2/) dzz^/^-^x +1-2Z + y)-'-"/^ 

dzz^l^-^x + 1-2Z + y)-'^/2| {X-1 + yfl\ 

Evaluating numerically the right hand side, as y ^ the curve F(x,y) approaches (Fig. 3) the limit function 
F{x,0) rising from F(1,0) = Tc to = lim^^oo i^(a;, 0) = Tc/X^, where X^o = {d - 2)/d [35]. Similarly, 
carrying out numerically the inversion of Eq. (70), the function x = f^^{t'^C, y) shows (Fig. 4) the approach to 
the singular behavior (41) in the limit y ^ 0. The parametric plot in Fig. 5 displays the approach of T{C,tw) 
toward T(C) of the form of Eq. (42) as tyj oo, which can be very slow if b is small. With the purpose of 
comparing further down with data for the XY model, in Fig. 6 we have also shown the parametric plot of the 
ZFC susceptibility. This has been obtained by plotting x(t,t^) against C{t,ty^,) for fixed f„ in two different 
cases: in panel (a), with d = 2.5 corresponding to & = 0.25, the approach to the asymptotic linear plot (45) is 
evident, while in panel (b), with d = 2.06 corresponding to the much smaller value b = 0.03, there is a much 
longer preasymptotic regime preceding the onset of the linear behavior. 

2. Quench below the critical line 

Using Eqs. (74) and (76) we may rewrite Eq. (35) as 

Teffit, t^, to) = Tx-'"^ + tt!^-^F{x, y){l - x-^/^) (79) 

where 

F{x,y) = M\Snr^^- . { fi^l±lV^''' . (80) 

^ ^ ' 4(a;'^/4-l) \x + l + yj ^ ^ 

In Fig. 7 we have plotted Tej / (t, t^, to) against the exact C{t, t^, to) obtaining T(C, t^), which shows the approach 
toward the form (49) of T(C) as t^ grows. 

3. Quench to T = with d = dL 

At d = 2 and T = 0, the second term in Eq. (55) vanishes and 60- — 0, yielding 

C{x,y) = 2x^/\x + l+y)-^ (81) 

while from Eqs. (79) and (80) 

Teff{t,t^,to) = (47r) (^|^i±^y . (82) 
Letting y ^ and eliminating x between the above two equations we obtain the explicit non trivial expression 

plotted in Fig.2c. We recall that the specific form of the function T(C) in the quench to {dL,T = 0) is 
nonuniversal. 
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VII. XY MODEL 



In the XY model the phase diagram of Fig.l includes the KT line and, as anticipated in the Introduction, in the 

d = 2 case the behavior of T(C) for < T < Tkt is expected to be given by Eq. (42). This is hard to detect from 
existing data, since the small values of 6 = r]{T)/z at the temperatures used in [39,40] would require to reach huge 
values of in order to observe the approach to (42). In Ref. [39] data are presented for the ZFC susceptibility. 
According to Eq. (44) asymptotically one should see the approach to the linear parametric plot (45) , as illustrated in 
Fig. 6a for the large- model. The data in Fig. 10 of Ref. [39] seem to be far from this behavior and, indeed, in [39] 
these data have been regarded as reminiscent of the non trivial parametric plot in the d = 3 Edwards-Anderson 
model. However, we believe that the non trivial pattern displayed by these data must be attributed to the large 
preasymptotic effect due to the small value of 6 = rj{T)/z = 0.03. Indeed, the pattern displayed by x(C, t^) in Fig. 10 
of Ref. [39] is strikingly similar to the one in Fig. 6b for the critical quench in the large- A' model, despite the fact that 
in the large- A^ model there is no KT transition. The data in Fig. 6b have been generated in order to reproduce the 
same value of h in the quench of the XY model considered in [39] . 

It would be interesting to have data for the XY model quenched at d^ and T = 0, where, according to the general 
argument of section III, x(C) is expected to display the non trivial behavior discussed in section IV. 

In the d = Z case in Ref. [41], the parametric plots of the FDR 

X{C)= lim ^ = 7F^ (84) 

t^^oo T,f^{c,t^) no 

are presented in the quenches to and below Tc- The data display the approach to the asymtotic parametric forms 

Y(r\- I ^°-=^c/T^ for C = 
^y^>-\ 1 for 0<C<1 



in the quench to Tc and 



, . _ , for C<M^ 
^^^> - 1 1 for M2 < C < 1 ^^^> 



in the quench below Tc, which correspond to Eqs. (42) and (49) for T{C). 

Furthermore, in the d = 3 model of particular interest is the measurement of the exponent a below Tc- The value 
a = 1/2 has been obtained both from the ZFC susceptibility [17] and from the direct mesurement of R{t,tw,to) [41]. 
This value gives additional and clearcut evidence in favour Eq. (54), which, in fact, predicts a = 1/2 when d = 3 and 
the order parameter is non conserved (z — 2) and vectorial {n — 2, du — A). Conversely, the conjecture [32] relating a 
to the density of defects, would predict a = 1, since in the XY model, with a vectorial order parameter, the density 
of defects goes like [10,42]. 



VIII. CONCLUSIONS 



In this paper we have overviewed the behavior of the effective temperature in the slow relaxation processes arising 
when systems, like a ferromagnet, with a simple pattern of ergodicity breaking in the low temperature state, are 
quenched from high temperature to or below Tc- On the basis of very general assumptions on the scaling properties 
of the autocorrelation and response functions, we have derived the scaling behavior of the effective temperature, 
which allows to account in a unified way for the different patterns displayed by T(C) in the different regions of the 
phase diagram. The primary results are i) in the critical quenches with Tc > 0, as a consequence of a = 6 > 0, 
T(C) displays always the form (42), which is trivial in the sense that T(C) = Tc except for the limiting value Too of 
Godreche and Luck at C = 0. In particular, this implies that also in the d = 2 XY model qucnc;lic(l to < T < Tkt 
T{C), and the equivalent parametric representations of the FDR or the ZFC susceptibility, must be trivial. The non 
trivial behavior reported in ref. [39] , then, must be regarded as preasymptotic. ii) The non trivial form of T{C) in the 
quenches to {d^^T = 0) requires a = 0. Since, due to weak ergodicity breaking, this process cannot be regarded as 
the continuation to Tc = of the critical quenches, a = cannot be explained on the basis of Eq. (11). Rather, the 
exponent a in the quenches below Tc must be expected to have a dimensionality dependence such that lim^^d^ a = 0. 
The phenomenological formula (54) fits quite well within this framework, although a deeper theoretical understanding 
of the problem is needed in order to put it on firmer grounds. 
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CI C 

FIG. 2. Parametric plot of the effective temperature T{C): (a) quench to Tc > 0, (b) quench to T < Tc (c) quench to 
{dL,T = 0). Panel (c) has been generated plotting the explicit form of T{C) given by Eq. (83) in the large- A'^ model. 
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FIG. 4. Approach to the limiting function x{C) of Eq. (41) for ty, = 10, 10^, 10^, 10* (top to bottom) in the critical quench 

of the 3d large- iV model. 
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FIG. 5. Approach of fe//(C,t„) toward T(C) of the form of Eq. (42) for tn, = 10, 10^, 10^, 10* (top to bottom) in the critical 
quench of the 3d large- model. 
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FIG. 7. Approach of feff{C,tyj) toward T(C) of the form of Eq. (49) for t-uj = 10^ 10^ 10^, lO'^ (bottom to top) in the 
quench to T = Tc/2 of the 3d large-TV model. 
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